© The Spatial Standard Observer 

Degrees of visibility and discriminability of targets in images can be estimated. 

Ames Research Center, Moffett Field, California 



The SSO Was Applied to These Two Images of targets viewed from a UAV. The left image contains a 
target calculated to have a visibility of 4.25 JND; the right image contains a similar target that has 
been reduced in contrast to have a visibility of 1 


The spatial standard observer is a com- 
putational model that provides a measure 
of the visibility of a target in a uniform 
background image or of the visual dis- 
criminability of two images. Standard ob- 
servers have long been used in science 
and industry to quantify the discrim- 
inability of colors. Color standard ob- 
servers address the spectral characteris- 
tics of visual stimuli, while the spatial 
standard observer (SSO) , as its name in- 
dicates, addresses spatial characteristics. 

The SSO is based on a model of human 
vision. The SSO was developed in a 
process that included evaluation of a num- 
ber of earlier mathematical models that 
address optical, physiological, and psy- 
chophysical aspects of spatial characteris- 
tics of human visual perception. Elements 
of the prior models are incorporated into 
the SSO, which is formulated as a compro- 
mise between accuracy and simplicity. The 
SSO operates on a digitized monochrome 
still image or on a pair of such images. The 
SSO consists of three submodels that oper- 
ate sequentially on the input image (s): 

1 . A contrast model, which converts an 
input monochrome image to a lumi- 
nance contrast image, wherein lumi- 
nance values are expressed as excur- 
sions from, and normalized to, a mean; 

2. A contrast-sensitivity-filter model that 
includes an oblique-effect filter (which 
accounts for the decline in contrast sen- 
sitivity at oblique viewing angles) ; and 

3. A spatial summation model, in which 
responses are spatially pooled by rais- 
ing each pixel to the power beta, 
adding the results, and raising the sum 
to the 1/(3 power. In this model, (3=2.9 


was found to be a suitable value. 

The net effect of the SSO is to compute 
a numerical measure of the perceptual 
strength of the single image, or of the vis- 
ible difference (denoted the perceptual 
distance) between two images. The unit 
of a measure used in the SSO is the just 
noticeable difference (JND), which is a 
standard measure of perceptual discrim- 
inability. A target that is just visible has a 
measure of 1 JND. 

The SSO was devised to satisfy an in- 
creasing need for a rapid, objective means 
of estimating degrees of visibility and dis- 
criminability of visual elements in scenes 
observed, not only by humans, but also by 
robotic vision systems, under a variety of 
circumstances. Examples of potential ap- 
plications of the SSO include evaluating 
vision from unpiloted aerial vehicles 
(UAVs) [see figure]; predicting visibility 


JND. 

of UAVs from other aircraft; estimating 
visibility, from a control tower, of aircraft 
on runways; measuring visibility, from a 
distance, of damage on aircraft and on a 
space shutde; evaluation of legibility of 
text, icons, or other symbols; specification 
of resolution of a camera or a display de- 
vice; inspection of display devices during 
manufacturing; estimating the quality of 
compressed digital video imagery; and 
predicting the outcomes of corrective 
laser eye surgery. 

This work was done by Andrew B. Watson 
and Albert J. Ahumada, Jr., of Ames Re- 
search Center. 

This invention is owned by NASA and a 
patent application has been filed. Inquiries 
concerning rights for the commercial use of this 
invention should be addressed to the Ames 
Technology Partnerships Division at (650) 
604-295 4. Refer to ARC-14569-1. 


Less-Complex Method of Classifying MPSK 

Nearly optimal performance can be obtained with less computation. 

NASA’s Jet Propulsion Laboratory, Pasadena, California 


An alternative to an optimal method of 
automated classification of signals modu- 
lated with M-ary phase-shift-keying (M-ary 
PSK or MPSK) has been derived. The alter- 
native method is approximate, but it offers 
nearly optimal performance and entails 
much less complexity, which translates to 
much less computation time. 

Modulation classification is becoming 


increasingly important in radio-communi- 
cation systems that utilize multiple data 
modulation schemes and include software- 
defined or software-controlled receivers. 
Such a receiver may “know” little a priori 
about an incoming signal but may be re- 
quired to correctly classify its data rate, 
modulation type, and forward error-cor- 
rection code before properly configuring 


itself to acquire and track the symbol tim- 
ing, carrier frequency, and phase, and ulti- 
mately produce decoded bits. Modulation 
classification has long been an important 
component of military interception of ini- 
tially unknown radio signals transmitted by 
adversaries. Modulation classification may 
also be useful for enabling cellular tele- 
phones to automatically recognize differ- 
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Performances of the ML and Approximate Methods were computed in simulations for a case of non- 
coherent reception of signals that had equal probability of being binary PSK or quaternary PSK. 


ent signal types and configure themselves 
accordingly. 

The concept of modulation classifica- 
tion as outlined in the preceding para- 
graph is quite general. However, at the 
present early stage of development, and 
for the purpose of describing the present 
alternative method, the term “modula- 
tion classification” or simply “classifica- 
tion” signifies, more specifically, a distinc- 
tion between M- ary and M'- ary PSK, 
where M and M' represent two different 
integer multiples of 2. 

Both the prior optimal method and the 
present alternative method require the ac- 
quisition of magnitude and phase values of 
a number (N) of consecutive baseband 
samples of the incoming signal + noise. The 
prior optimal method is based on a maxi- 
mum-likelihood (ML) classification rule 
that requires a calculation of likelihood 
functions for the M and M' hypotheses: 
Each likelihood function is an integral, over 
a full cycle of carrier phase, of a complicated 
sum of functions of the baseband sample 
values, the carrier phase, the carrier-signal 
and noise magnitudes, and M or M\ Then 
the likelihood ratio, defined as the ratio be- 
tween the likelihood functions, is com- 
puted, leading to the choice of whichever 
hypothesis — M or M 1 — is more likely. 

tn the alternative method, the integral 
in each likelihood function is approxi- 
mated by a sum over values of the inte- 


grand sampled at a number, 7, of equally 
spaced values of carrier phase. Used in 
this way, / is a parameter that can be ad- 
justed to trade computational complex- 
ity against the probability of misclassifi- 
cation. tn the limit as / — > oo, one 
obtains the integral form of the likeli- 
hood function and thus recovers the ML 
classification. 

The present approximate method has 
been tested in comparison with the ML 


method by means of computational simu- 
lations. The results of the simulations 
have shown that the performance (as 
quantified by probability of misclassifica- 
tion) of the approximate method is nearly 
indistinguishable from that of the ML 
method (see figure). 

This work was done by Jon Hamkins of Cal- 
tech for NASA’ s Jet Propulsion Laboratory. 
Further information is contained in a TSP (see 
page 1 ). NPO-40965 


Improvement in Recursive Hierarchical Segmentation of Data 

Goddard Space Flight Center, Greenbelt, Maryland 


A further modification has been made 
in the algorithm and implementing soft- 
ware reported in “Modified Recursive Hi- 
erarchical Segmentation of Data” (GSC- 
14681-1), NASA Tech Briefs, Vol. 30, No. 6 
(June 2006), page 51. That software per- 
forms recursive hierarchical segmentation 


of data having spatial characteristics (e.g., 
spectral-image data) . The output of a prior 
version of the software contained artifacts, 
including spurious segmentation-image re- 
gions bounded by processing-window 
edges. The modification for suppressing 
the artifacts, mentioned in the cited article, 



A Segmentation of the Landsat ETM+ Image displayed on the left is shown on the right. The new ap- 
proach eliminates processing artifacts. 


was addition of a subroutine that analyzes 
data in the vicinities of seams to find pairs 
of regions that tend to lie adjacent to each 
other on opposite sides of the seams. 
Within each such pair, pixels in one region 
that are more similar to pixels in the other 
region are reassigned to the other region. 
The present modification provides for a 
parameter ranging from 0 to 1 for control- 
ling the relative priority of merges between 
spatially adjacent and spatially non-adja- 
cent regions. At 1, spatially-adjacent-/spa- 
tially-non-adjacent-region merges have 
equal priority. At 0, only spatially-adjacent- 
region merges (no spectral clustering) are 
allowed. Between 0 and 1, spatially-adja- 
cent-region merges have priority over spa- 
tially-non-adjacent ones. 

This program was written by James C. 
Tilton of Goddard Space Flight Center. 
Further information is contained in a TSP 
(see page 1 ). 

GSC- 14 994-1 


NASA Tech Briefs, September 2006 


51 


